Anomalous diffusion at the Anderson transitions 



Tomi Ohtsuki 

Department of Physics, Sophia University, Kioi-cho 7-1, Chiyoda-ku, Tokyo 102, Japan 

Tohru Kawarabayashi 

Institute for Solid State Physics, University of Tokyo, Roppongi, Minato-ku, Tokyo 106, Japan 

(February 1, 2008) 

Abstract 

Diffusion of electrons in three dimensional disordered systems is investigated 
numerically for all the three universality classes, namely orthogonal, unitary 
and symplectic ensembles. The second moment of the wave packet < r 2 (t) > 
at the Anderson transition is shown to behave as ~ t a {a ~ 2/3). From the 
temporal autocorrelation function C(t), the fractal dimension D2 is deduced, 
which is almost half the value of space dimension for all the universality 
classes. 

Metal-insulator transitions are one of the most extensively investigated subjects in con- 
densed matter physics. Especially interesting is the quantum phase transition, where the 
transition is driven by changing the parameter of quantum systems instead of temperature. 
The Anderson transition |],|| is a typical example, where extended electronic states become 
localized with the increase of disorder. 

Much effort has been devoted to clarify the Anderson transition, both experimentally and 
theoretically. In the metallic regime where the electronic states are extended, the transition 
is defined by the vanishing conductivity a as the strength of disorder W is increased. It is 
characterized by the critical exponent s as 

o-~(W c -W)\ (1) 



with W c the critical disorder. In the insulating regime where the states are localized, it is 
most clearly seen by the divergence of the localization length £i oc as 

6oc ~ (w - w c yr (2) 



From the one parameter scaling theory, the exponent v is related to s by the Wegner's 
scaling law, || 

s = {d- 2)v, (3) 

d being the dimensionality of the system, and once v is determined, we can predict the 
behavior of the conductivity near the Anderson transition. 

The behaviors of the localization length and the conductivity are conjectured to be 
universal, i.e., v does not depend on the detail of the system. They are determined only by 
the basic symmetry of the system under the operation of time reversal or spin rotation. |||6[] 
Systems with time reversal and spin rotation symmetry are called orthogonal ensemble, while 
those with only time reversal symmetry are called symplectic ensemble. Systems without 
time reversal symmetry are unitary ensemble. 

The value of v has been calculated for three dimensional system by using the finite-size 
scaling argument. It is estimated to be 1.4 ±0.1 for orthogonal ensemble, [[| 1.3 ± 0.2 for 
unitary ensemble, [§-[K| an d again 1.3 ± 0.2 for symplectic ensemble. |Ll,r^| These facts 



indicate that the critical behavior of conductivity as well as the localization length does not 
depend significantly on the symmetry of the system. 

On the other hand, recent analyses on the energy level statistics at the Anderson transi- 
tion show that the level statistics do depend on the symmetry of the system, though they are 
independent of system size or model. |TT|-F2Tf This universal behavior is related to the scale 
invariance at the transition, where eigenfunctions show fractal structure. Peculiar behavior 
just at the transition is now attracting a lot of attention. 

In this paper, we numerically discuss electron diffusion at the Anderson transition in 
three-dimensional (3D) disordered systems. The diffusion coefficient becomes size dependent 



at the transition, which leads to the increase of the mean square diffusion length < r 2 (t) > 



as t a with the exponent a < 1. |22|] The return probability also decreases as power law, 
reflecting the fractal structure of the wave function. |23| , |24|| 

First we discuss the behavior of < r 2 (t) >. As in the percolation theory, p5[ let us 
assume the scaling form 

< r 2 {t) >= Ct kl f {{W c - W)t k2 ) . (4) 

In the metallic regime, < r 2 (t) > increases as 2dDt where D is the diffusion constant. It is 
related to the conductivity from the Einstein relation, and behaves as (W c — W) s . Therefore, 
f(x) should be proportional to x s in the limit of large x. 

In the insulating regime, the wave packet ceases to diffuse if the diffusion radius becomes 
the order of localization length. Therefore we have < r 2 (t) >~ £f oc ~ (W c — W)~ 2u and 
f(x) is proportional to (—x)~ 2u when —x is large enough. From this argument, we have 

ki + sk 2 = 1 , 

(5) 

k\ — Ivki = , 



and consequently 



k\ = , k 2 = . (6) 



Using the scaling relation s = (d — 2)z/, we have 



We therefore expect 



<r 2 {t)>~t 2 '\ (8) 



at the Anderson transition {W = W c ). 

At the transition, the wave function has fractal structure. In this situation, if the radius 
of the wave packet is r, the return probability P(t) is proportional to r~ D ' 2 . From (^|), we 
have 



P(t) ~ r D2 / d . (9) 

From the behavior of P(t), we can determine the fractal dimension D 2 . This intuitive 
argument agrees with the more detailed analysis of the scaling behavior of the dynamical 



diffusion coefficient D(q,u). f23| , |24| 



In the actual simulation, we have adopted tight binding Hamiltonian 

H = £ Vi^jycl^cj^ + Wic{ a c i)a , (10) 

<i,j>,a,a' i,<J 

where i,j denote the lattice site, and a, a' the spin. In the orthogonal case, V ita .j i(T i = V5 a>a i is 
real,while Vi^j^i is V exp(i(pi y j)5 cr y with </>y the Peierls phase factor in the unitary case. In 
both cases, no spin flip process is included. In the symplectic case, the hopping is described 
by 

Vi^i-ky = V[exp(-i9a k )] a y, k = x,y, z, (11) 



where oVs are Pauli matrices. [0,0] We have assumed the simple cubic structure with the 
lattice constant taken to be unity. Only the nearest neighbor coupling is assumed. The site- 
diagonal potentials W,- t are assumed to be distributed independently, and their distribution 
is taken to be uniform in the range [-W/2, W/2]. 

Instead of diagonalizing the system directly, we solve numerically the time- dependent 
Schrodinger equations. |27| , |28|| We evaluate the time-evolution operator U(t) = exp(—iHt/h) 



by using the decomposition formula for exponential operators. p9[ The n-th order decom- 
position U n satisfies the condition 

U(8t) = U n {5t) + 0{5t n+1 ). (12) 

We have adopted the same forth-order decomposition, as in the previous papers, p7|j28 



given by 

u 4 = u 2 (-i P t/h)u 2 (-i(i - 2 P )t/h)u 2 (-ipt/h) (13) 
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with 

JJ 2 (x) = e xHl ^ 2 ■ ■ ■ e xH <i-i/ 2 e xH <i e xH <i-i/ 2 . . . e x ' Hl / 2 , 

where H = Hi + ■ • ■ + H q and p = (2 — v^) -1 . The decomposition is made so that each 
Hamiltonian Hi(i = 1, ■ • • , q) should consist of commuting terms. 

The actual simulations have been done in systems with 59 x 59 x 59 lattice sites for 
orthogonal and unitary ensembles, while 69 x 69 x 69 lattice sites are used to discuss the 
symplectic case. In each case, average over 10 independent realizations of random potentials 
has been performed. The initial wave packet is build by diagonalizing a spherical system 
with radius R — 3 located at the center of the whole system. We use the wave packet whose 
energy is closest to the band center. 

To discuss the properties at the vicinity of the transition, we have set W = W c = 16.5V 
for orthogonal case. || In the unitary case we assume that the magnetic field is parallel 
to the ^-direction, and the magnetic flux penetrating the x-y plane unit cell is set to be 
0.1 times the flux quantum. The resulting critical disorder is W = W c = 17.8V. || For 
the symplectic ensemble, we have set 9 = tt/6 in ([11]) and W is again set to the critical 
value W c = 19.0V. JLT| The time step is chosen to be St = 0.2/i/Vi where Vi is the hopping 
amplitude without spin flip process. 

We first discuss the critical behavior of the second moment of the wave packet < r 2 (t) > 
defined as 

< r 2 (t) > c =< t\r 2 \t > - < t\x\t > 2 - < t\y\t > 2 - < t\z\t > 2 , (14) 

where \t > denotes the state at time t. In Fig. 1, we plot it as the function of time tV/h for 
the three universality classes. The solid line corresponds to the orthogonal case, the broken 
line to the unitary, and the dotted line to the symplectic one. The standard deviations with 
respect to 10 realizations of random potential configurations are typically less than 5%. From 
these behavior, < r 2 (t) > c is estimated to increase as t a , with a = 0.67±0.02 for orthogonal 
case, 0.66±0.02 for unitary case and 0.69±0.02 for symplectic case. The excellent agreement 
of a with 2/d = 2/3 confirms the scaling form (£|) and the scaling relation s = (d — 2)u. 
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[ Fig. 1 ] 

Now we discuss the return probability. In Fig. 2, we plot the temporal autocorrelation 
function []3D| C{t) defined by the overlap function between the initial state and the state at 
time t as 

C(t) = - f dt'| < t'\0 > | 2 = - f dt f P(t'). (15) 
t Jo t Jo 

Average of logC(t) over 10 random potential configurations has been performed, and the 

standard deviation is indicated for orthogonal case (the standard deviations for the other 

cases are not shown but they are almost the same). From the asymptotic behavior, we 

estimate the fractal dimension D 2 as 1.5 ± 0.2 for orthogonal, 1.7 ± 0.2 for unitary and 

1.6 ± 0.2 for symplectic case. In the case of orthogonal ensemble, the value 1.5 ± 0.2 agrees 

with the results obtained previously by the direct diagonalization. P^pT|J52] ] Our new results 

for unitary and symplectic cases show that the fractal dimensionality D 2 does not depend 

strongly on the symmetry. The results are summarized in Table I. 

[ Fig. 2] 

In conclusion, we have studied the diffusion of electron in 3D disordered systems at 
the Anderson transition by numerically solving the Schrodinger equation. The anomalous 
diffusion r 2 (t) ~ t 2 / 3 has been clearly observed, which is expected from the scaling form 
(|j) and the scaling relation s = (d — 2)u. The fractal dimensionality D 2 is also estimated. 
For all three universality classes, D 2 is almost half the space dimension. It is interesting 
to note that the values D 2 for two dimensional Anderson transitions, namely the quantum 



Hall and symplectic systems, are almost the same as well (1.62 ± 0.02 for the former p3|j34 



and 1.66 ± 0.05 for the latter. |28| . |3"5|j ). As discussed by Brandes et al, |24| , |3"6|| the fractal 
properties are reflected in the temperature dependence of the inelastic scattering time r in 
at the Anderson transition in relatively high temperature. Our results for unitary and 
symplectic cases indicate that almost the same temperature dependencies are observed in 
all 3D Anderson transitions. 
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Finally, let us discuss the temperature dependence of the conductivity at the transition. 
By the anomalous diffusion, the relation between the inelastic scattering length Z in and the 
inelastic scattering time is modified to be 

l in oc if/ 3 . (16) 

Then the effective diffusion constant D e R observed at finite temperatures is 

D cS ~ 4 „ T r Q y\ (17) 

leading to a ~ Tin . At sufficiently low temperature, r-' 1 is proportional to the temperature 
T. Experimentally observed T 1//3 behavior of the conductivity at the transition is thus 
consistent with the present scaling argument using the scaling relation s = (d — 2)v. It is 



recently suggested |37] that the scaling relation is modified in the interacting system where 
Anderson-Mott transition occurs. Careful investigation of the temperature dependence of 
the conductivity at the transition will clarify the nature of the transition. 

The authors are grateful to Professor Yoshiyuki Ono and Dr. T. Brandes for fruitful 
discussions. This work is in part financed by the Grants-in-Aid 08740327 from the Ministry 
of Education, Science and Culture. The numerical calculations have been in part performed 
on a FACOM VPP500 of Institute for Solid State Physics, University of Tokyo. 
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TABLES 



TABLE I. 

classes. 


Summary of the exponent 


a and the fractal dimensionality D 2 


for three universality 




orthogonal 


unitary 


symplectic 


a 


0.67 ±0.2 


0.66 ±0.2 


0.69 ±0.2 


D 2 


1.5 ±0.2 


1.7 ±0.2 


1.6 ±0.2 



Figure captions 

Fig. 1: The growth of the second moment < r 2 (t) > c of the wave packet. The solid line 
corresponds to the orthogonal case, the broken line to the unitary and the dotted line to the 
symplectic one. In large t regime, t 2//3 behavior is clearly seen. 

Fig. 2: The time- dependence of the auto-correlation function C(t) . The bars around 
the data for orthogonal case indicate the standard deviation with respect to 10 realizations 
of random potential configurations. They are almost the same for all universality classes, so 
only those for the orthogonal case are shown for simplicity. 
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Fig. 1 Ohtsuki and Kawarabayashi 




Fig. 2 Ohtsuki and Kawarabayashi 




